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Abstract 

The problem of the hierarchy of lifetimes of charmed hadrons is 
reviewed. The QCD-based theory of preasymptotic effects in inclu- 
sive weak decays dating back to the beginning of the eighties is now 
entering its mature phase. Combining recent and old results we ar- 
gue that the observed hierarchy reflects most intimate features of the 
hadronic structure. The problem of a wide spread of lifetimes of 
charmed hadrons is addressed. We speculate on what is to be expected 
from QCD to provide the observed pattern. A number of predictions 
is given for the hierarchy of lifetimes in the family of the beautiful 
hadrons. 



1 Introduction 



From early days of QCD it is well-known that the lifetimes of all particles Hq 
containing a heavy quark Q and light quarks are the same and equal to that 
of a 'free' heavy quark in the limit itlq — >■ oo {mq is the mass of the quark Q). 
At the same time the experimentally observed total widths in the charmed 
family differ by as much as an order of magnitude for and This 
scatter clearly shows that mq ~ 1.4 GeV is far below the asymptotic domain, 
and the questions arise whether the preasymptotic effects are understood in 
the present-day QCD and what one should expect for the h quark family. 
These questions are addressed in the present talk. As we will see later, the 
hierarchy of the lifetimes presents a junction of different aspects of the theory 
of strong interactions. The two main elements are (i) a systematic expansion 
of amplitudes in the inverse mass of the heavy quark; (ii) estimates of matrix 
elements of local operators (appearing in the expansion) over the hadronic 
states Hq. The first step, the construction of the rriq^ expansion, takes into 
account the quark-gluon dynamics at short distances and is clean in the 
sense that it is based on fundamental QCD and is, thus, model-independent 
(with a certain reservation, to be discussed below). The technology of m'q^ 
expansions experienced a dramatic development in recent years |T|, ^. The 
corresponding achievements are incorporated in our discussion. As for the 
hadronic matrix elements, they reflect the hadron structure at large distances. 
As usually, here the progress was limited, if at all, and our purpose is to 
outline some possible approaches and reveal the impact of this or that specific 
mechanism on the lifetime hierarchy. 

The 'scientific' approaches to the issue of the lifetime hierarchy date back 
to the beginning of the eighties. It was pointed out that the total widths 
are related, by means of unitarity, to the imaginary part of certain forward 
scattering amplitudes |0] . This simple and rather obvious observation paved 
the way to a generalization of the Wilson operator product expansion (OPE) 
allowing one to represent all amplitudes of interest in the form of a series 
in rrtg^ [H, |, ^, an expansion essentially equivalent to what is called now 
the Heavy Quark Effective Theory (HQET) 0, |^ (for reviews see e.g. ref. 
[0]). Shortly after it was discovered that the terms behaving like powers of 



niq should be supplemented by hybrid logarithms |]T0|- These logarithms 
are a routine part of the present-day HQET. The analysis of preasymptotic 
effects in charm and beauty combining both power series in m'^ and the 
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hybrid logarithms has been first carried out in in the mid-eighties. At 
those days the data on baryons were scarce even in the c quark family, and 
next to nothing was known about the lifetimes in the b family. It is not 
surprising, therefore, that the emphasis in these works was on the operators 
which differentiate between the and D^. These operators necessarily 
involve the spectator light quark and have dimension 6 or higher. A set 
of dimension 6 operators was analyzed in P] Q with definite predictions 
concerning the hierarchy of the lifetimes. Now, ten years later, it became 
clear that the analysis has to be updated for many reasons. 

First, the numerical values of some crucial parameters turned out to be 
higher than previously thought. As a result, a typical scale of the preasymp- 
totic 'corrections' becomes larger to the extent that the 'corrections' over- 
shadow the leading term in the charmed family (so that, strictly speaking, 
they can not be called corrections anymore). Even the beauty family is still 
not quite in the asymptotic domain. 

The enhanced role of the preasymptotic effects leads, in turn, to the 
necessity of including other operators, not fully considered in this problem 
previously. There are two operators of dimension 5 and one new operator 
of dimension 6. Understanding the ratio of the baryon-to-meson lifetimes is 
impossible without proper consideration of these operators. 

As we will see below, for rric ~ 1.4 GeV the convergence of the expansion 
is unexpectedly so bad that even operators of dimension 7 and higher must 
play an important role. Not much can be said about them at the moment, 
which means, of course, that all results referring to the c quark family are 
doomed to be valid at best qualitatively, if at all. We have gloomy suspicions, 
to be shared with everybody, that the QCD-based 'inclusive' approach will 
never become fully quantitative in the charmed family. In the b quark family 
the situation seems to be more favorable. Here the the expansion parameter 
m^^ is sufficiently small to guarantee a reasonable convergence of the uiq^ 
series, and limiting to operators of dimension 5 and 6 is quite reasonable. 

Next, one can use symmetry arguments and models to estimate the 
hadronic matrix elements of these operators, establishing in this way a hier- 
archy of lifetimes. 

This stage, evaluation of the matrix elements, is the most vulnerable point 
of the existing theory, its Achilles' heel. One must frankly admit that at least 

^Ref. ignores the existence of the hybrid logarithms. 
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some of the relevant matrix elements are not controlled theoretically today. 
Hence, one is forced to move onto swampy soil of models. We invoke a prim- 
itive quark model for orientation. What is the outcome? Honestly, it should 
be admitted that the charmed family is too light for serious theorists building 
1 / tuq expansions. With some imagination, however, we still dare to say that 
the pattern of the charm lifetimes is understandable. And further efforts, 
both experimental and theoretical, are worthwhile - the observed hierarchy 
(once it is actually established) would imply certain relations between the 
hadronic matrix elements of the operators involved, a unique information, 
inaccessible in any other way. 

Using the fact that this talk is not subject to the scrutiny of refereeing we 
allow ourselves to indulge in a fantasy treating the l/rric expansion and the 
primitive quark model mentioned in a rather speculative way. An element of 
wishful thinking is quite obvious in this talk, and we bring our apologies for 
this fact. Our only justification is a hope that the corresponding speculations 
will catalyze a more 'scientific' analysis in future. 

We start with a necessary theoretical background and a brief historic 
review. Then we proceed to the picture of the c quark family. Finally, a few 
remarks concerning our expectations in the b quark family are given. 



2 Theoretical Excursion 
2.1 General picture 

In prehistoric times (the seventies) the inclusive heavy flavor decays were 
treated in the most straightforward way - the width of the Q containing 
hadron Hq was assumed to be equal to the probability of the free Q quark 
decay. An immediate consequence is the equality of the widths of all hadrons 
belonging to the Q family. 

Although in the limit rriQ oo this prescription is certainly correct a 
single glance on what we actually have in the charmed family (see Table 1) 
shows how far the existing quarks are from the asymptotic domain and how 
important preasymptotic corrections are. Even leaving aside S° where the 
data are rather poor and the error bars are large we can say that the lifetimes 
span the interval from roughly 10 x 10"^^ sec for (the most long-living 



charmed particle) to 2 x 10 for Ac (see ref. |TT| for the recent analysis of 
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the experimental data). There is a drastic imparity between and and 
between mesons and baryons. Do we understand this hierarchy theoretically? 

The idea lying at the base of the modern theory of the preasymptotic 
effects is simple and can be compared to a well-known text-book problem. 
We mean the problem of the nuclear (3 decay. If the energy released in the (3 
electron is much larger than a typical binding energy of both electrons in the 
shell and nuclcons in the nucleus (the latter is unrealistic, of course) one can 
forget about the electrons in the shell and other bound state effects in calcu- 
lating the decay probability. In the high precision calculations, however, the 
bound state effects should be accounted for. In particular, the soft electrons 
in the shell do interfere with the j3 electron (the Pauli exclusion principle), 
and this interaction generates corrections proportional to inverse powers of 
the energy release. 

2.2 The inverse quark mass expansion 

Likewise, in the heavy quark decay all participating degrees of freedom can 
be either hard ( the corresponding energy is of order mg) or soft. In the 
leading approximation we ignore the soft degrees of freedom and treat the 
hard ones as free - a good old free quark decay (Fig. 1). The inclusion 
of the soft degrees of freedom in analysis generates power non-perturbative 
corrections. An obvious example of the soft degree of freedom is the spectator 
light quark. This is not the end of the story, of course. Soft gluons in the 
initial light cloud surrounding Q or those emitted by the final light quarks is 
another example (see Fig. 2). 

The formal basis allowing one to systematically catalog the preasymptotic 
non-perturbative corrections is provided by Wilsonian operator product ex- 
pansion (OPE). More exactly, what is readily calculable are the so called 
condensate non-perturbative effects. 

Although the word OPE might sound frightening to non-experts, con- 
ceptually this is a very simple formalism, just a procedure of a systematic 
separation of large and short distance contributions in QCD (hard and soft 
degrees of freedom). We further assume that the short-distance contribu- 
tion residing in the coefficient functions can be found perturbatively, while 
all non-perturbative effects are attributed to matrix elements of the oper- 
ators appearing in the expansion. This procedure is sometimes called the 
practical version of OPE. We will refer to it as to the 'standard'. It is clear 
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that this standard prescription is by no means exact since in principle there 
are non-perturbative contributions even at short distances modifying the 
coefficient functions. (They will be referred to as hard or non-condensate 
non-perturbative terms). The latter are poorly controllable theoretically 
at present, and we will just assume that they can be disregarded for the 
time being. It seems quite probable, though, that they are non-negligible 
in the charmed family. Anyhow, we will disregard them. The 'hard' non- 
perturbative terms should be (and probably will be) targeted in future anal- 
yses of the preasymptotic effects in the charm decays. 

Constructing the operator product expansion we obtain the inclusive de- 
cay probabilities as an expansion in m'^ . Practically all terms up to m'^ are 
known at present. We will discuss them and then speculate on higher order 
terms. 

The statement above is an oversimplification. Apart from the power de- 
pendence on m'^ one encounters logarithmic dependencies on this parameter 
- the so called hybrid anomalous dimensions 0] (matching logarithms of 
HQET). For pedagogical purposes let us forget for a while about the hybrid 
logs. We will comment on them later. 

We start from the effective weak lagrangian describing non-leptonic de- 
cays 

£(/.) = ^VscVUciOi + C2O2) (1) 

where Ci 2 are Wilson coefficients accounting for the loop momenta from ix 
up to Mw and known in perturbation theory while Oi^2 are operators, 

0^ = {sT^c){uT^d), O2 = (s,r^c')(M,r^rf0- (2) 

Here = 7^(1 + 75). Eqs. (|l]), @) present the lagrangian relevant to 
the c — >■ sdu transition. The penguin graphs showing up at the 1% level 
are omitted. The coefficients Ci^2 are known in the leading log and next-to- 
leading approximations [|1^, [1^ ; 



ci = ^(C+ + C_), C2 = ^(C+-C_). (3) 



In the leading approximation 



d± 
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The numerical values of the coefficients depend on the choice of param- 
eters. The physical quantities, like the non-leptonic widths, can not depend 
on the arbitrary scale /z. The /z dependence explicit in the weak lagrangian 
(|l]) is compensated in the process of calculation of F's. In the leading log 
approximation we limit ourselves to one can set /i = from the very be- 
ginning. As for Aqcd the current fashion bends towards a rather high value. 



Aqcd ~ 300 MeV |Tj]. Our gut feeling is that the genuine value of Aqcd is 
actually smaller, but this is no time to indulge in the discussion of this issue. 
Accepting this estimate we get 

C+ ~ 0.74, ~ 1.8. (5) 

Now, we construct the transition operator T{c X ^ c), 

f = ij d^xT{C{x)Cm = J2 C^O, (6) 

i 

describing a diagonal amplitude with the heavy quark c in the initial and 
final state (with identical momenta). The transition operator T is built by 
means of OPE as an expansion in local operators Oi. The lowest-dimension 
operator in T(c —> X c) is cc, and the complete perturbative prediction 
corresponds to the perturbative calculation of the coefficient of this operator. 
In calculating the coefficient of cc we treat the light s, d, u quarks in X as 
hard and neglect the soft modes. Say, we ignore the fact that in a part of the 
phase space the u quark line is soft and can not be treated perturbatively. 
Likewise, we ignore interaction with the soft gluons. 

Our task - the task of the theory of preasymptotic corrections - is the 
analysis of the influence of the soft modes in the quark and gluon fields 
manifesting themselves as a series of high-dimension operators in T. 

Once the expansion is built we average T over the hadronic state of 
interest, charmed mesons and baryons in the case at hand. At this stage the 
non-perturbative large distance dynamics enters through the matrix elements 
of the operators of dimension 5 and higher. (There are no operators of 
dimension 4). 

Finally, the imaginary part of < Hc\T\Hc > presents the He width sought 

for, 

= ^< H,\lmf{c ^X^c)\H,>. (7) 



He 
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For experts we hasten to add a remark concerning peculiarities of the kine- 
matics of the amphtude under consideration. It is essentially Minkowskian. 
Still the operator product expansion originally formulated for deep Euclidean 
kinematics can be used. It is important that we take the full imaginary part 
corresponding to all possible cuts of the diagrams (totally inclusive final state 



X). Omitting some of the cuts would lead to infrared unstable results ^5 
Since the decay rates we calculate present integrated quantities over the in- 
terval of energies up to rric the corresponding prediction for F's, through 
analyticity, is related to the calculation of the OPE coefficients with the 
characteristic off-shellness ~ rric. For comparison let us mention a similar 
analysis of the r decays [0. 

2.3 Catalog of relevant operators 

The leading operator in the expansion of T is 

Oo = cc. (8) 

Keeping only cc is equivalent to perturbative calculation of the coefficients 
Ci^2 followed by a perturbative calculation of the decay width F through the 
relation 

F^^"* = — Im < elf |c>, (9) 
rric 

where taking a 'mathematical' matrix element over the c quark singles out 
the operator cc. 

It is worth noting that the contribution of Oq in the physical width F 
(not to be confused with r^^*"*!) is given by 

F = 2ImCo^^ < HJOolHc > . 

The matrix element {1/2MhJ < Hc\Oq\Hc > is 1 plus (9(m~^) corrections 
(see eq. (12) below). 

The next operator in our catalog is 

1 -* 
Og = -ca^^G^^c -caBc, (10) 

with dimension 5. This operator as well as O^^, see below, generates a differ- 
ence between mesons on one hand and baryons, on the other. 
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Let us make a side remark concerning another dimension 5 spin opera- 
tor, 

a = c[D' - {vDf]c, (11) 

where v is the four-velocity of the hadron He- This operator does not appear 
exphcitly in the calculation of the total widths due to the fact that 0-,^ is not 
Lorentz scalar. On the other, the relation [|l], ^ 

cc = -c^,v,c - + + -^f-ci^t^cY^qi^t^q + 0(l/m^) (12) 

used in the calculation of the matrix element < Hc\cc\Hc > contains 0-,^. 
(Here are the color generators and g is a generic notation for light quarks). 

We pass now to the 0{m~^) terms due to operators of dimension 6. On 
general grounds one can limit the number of these operators to a few following 
structures |T^. 

(i) Four- quark operators 
These have the generic form 

O^g = (crg)(grc) 

where g is a light quark field and F presents, in the case at hand, a combi- 
nation of the Lorentz and color matrices. The coefficients of the four-quark 
operators have been first calculated in refs. 0, ||, |]. They are determined 
by the one-loop graphs of the type presented on Fig. 3. 

The four-fermion operators are responsible for the splittings between 
and and the splittings between different charmed baryons. 

(ii) Dimension 6 quark-gluon operators 

Apart from the heavy quark fields c and c relevant operators contain 
the gluon field strength tensor G^^p (they may also include the covariant 
derivatives). These operators are generated by two-loop graphs (fig. 2) and, 
hence, their coefficients are numerically strongly suppressed compared to the 
four-quark case above. 

By using equations of motion it is not difficult to show [O] that there 



are only two options for the spin-zero quark-gluon operators of dimension 6, 
namely 

and 

Oe = ca^,yGf,p-fpiD,yC ^ ca E X iDc, (13) 
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where E is the chromoelectric field. All other operators of dimension 6 with 
the gluon field are reducible to the above. The operator c{D ^u)^ vC is 
actually a four-quark operator since 

Its coefficient contains extra ag/Tt., however, compared to the four-quark op- 
erators coming from the one-loop graphs, and the corresponding contribution 
can be neglected. 

As for Oe its contribution turns out to be small numerically compared 
to that of Og, so that it plays no role in what follows. We mention it for 
completeness and will essentially forget about Oe from now on. 

For pedagogical purposes it is instructive to write down the decomposition 
of the operators above in terms of the HQET field, 

h,{x) = e-~t^^c(x), (14) 

although, certainly, this decomposition adds nothing in practical terms. In 
eq. (|T^) c is the normal Dirac bispinor, and v is the four-velocity of the heavy 
hadron; v = •y^'Va- 

Then, say, for cc we essentially repeat eq. (12), 

Co hc'^^v^hc - ^~2^<^^'^^c ~ 2rn?^''^^^'' "' ^^^^ 

and 

Og ^ -hcaBhc + ... (16) 

It is possible to prove that all dimension 6 operators appearing in eqs. (|T5|), 
(|TB|) can be reduced to four-quark operators suppressed by They can be 
consistently discarded. 

2.4 Coefficients 

The operators listed above appear in T each with its own coefficient. Practi- 
cally all of them have been calculated previously - for Oq in the prehistoric 
times. The coefficients of the operators of dimension 5 were found in the 
recent works [|l], |^ and for the spectator-sensitive operators of dimension 6 
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in refs. We supplement this set by a new spectator-blind operator of 

dimension 6, and consider a new source of l/mq corrections - the expansion 
of matrix elements of operators of dimensions 3 and 5 in terms of l/mq. 
The prehistoric coefficient Cq is obviously equal to 

ImCo = i(3ro)r/ (17) 



where 

^° = 192vr3 ^^^^ 
is a very convenient unit for all widths to be discussed here. The factor r] 
reflects the hard gluon exchanges hidden in C±, 

7] = i « 1.5. (19) 

If the hard-gluon exchanges are accounted for in the next-to-leading order 



where the explicit expression for J is given in ref. |T8|, see also |T^. Since 



the picture of the lifetime hierarchy we are aimed at is qualitative at best we 
limit ourselves to the leading approximation. 

The strange quark mass is neglected, here and below, systematically. 
Since we deal exclusively with the current quarks (there is no such notion as 
the constituent quark in the QCD-based approach). 



ml /ml ~ 10^2 

and the approximation m^ = m^ = m^ = in calculating the inclusive widths 
seems to be not bad. One should be alert, however - possible deviations from 
this limit might be larger than expected, and an estimate of the effects due 
to non-vanishing M'^ is deflnitely in the short list of questions for further 
investigation. 

Now we are approaching a more interesting and less trivial stage of the 
analysis - we pass to soft modes. (A reminder: we discard non-condensate 
non-perturbative effects. This is not because they are necessarily small in 
the charm decays but because we do not know how to treat them. Thus, here 



10 



is another item from the short hst - hopefully, a courageous theorist can be 
found for its investigation.) The first effect to show up is that due to the 
chromomagnetic operator caBc describing the correlation of the heavy quark 
spin in He with the chromomagnetic field in the light cloud surrounding c. 
Its dimension is 5, hence the corresponding contribution in F is 0{m~'^). 

The chromomagnetic field B appears because both the initial quark c and 
the light fast quarks produced in the decay are coupled to it - they do not 
live in the empty space but, rather, in a medium, the 'brown muck' around 
the heavy quark (we prefer to call this medium 'light cloud'). There are two 
sources of Oq'- the direct coupling to d and the coupling to c. In the latter 
case Og appears through the equations of motion for the c field. We will not 



dwell on the issue referring to the original publications [|1], ^, |T^ for further 



details. To the leading order in as the coefficients Cq and Ct, are 

Cg = -To\{9,cI-2cV), (20) 

a = 0. (21) 

appears only when we take the matrix element of (9o, along with an 
additional piece with Cg, see eq. (12). 

Although both operators, Oq and (9,^, are of dimension 5 their matrix ele- 
ments over Hf. may contain also corrections 0{m~^) corresponding to 0{m~^) 
terms in F's which turn out to be insignificant numerically. 

We now move to the operators of dimension 6. As we argued above 
there are only two classes of operators that we have to consider here. All 
other operators can be reduced (through equations of motion) to color blind 
four-quark operators, whose Wilson coefficients are additionally suppressed 
by Us- Since considering corrections to the Wilson coefficients is beyond 
our accuracy now we shall neglect them. We are left with two classes of 
operators. First, and most important, are spectator-sensitive four-quark op- 
erators, whose analysis was carried out in refs. ^, |^ . Their coefficients are 
numerically enhanced since they are generated by one-loop graphs. These 
operators lead to the lifetime differences between different mesons and differ- 
ent baryons. There are three distinct physical mechanisms summarized by 
the four-quark operators: annihilation, interference and quark-quark scat- 
tering in the baryons. Technically the distinction manifests itself in the fact 
that each of the three light quark lines, d oi s can be soft. Pictorially 
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the softness of the hne is depicted as follows: instead of drawing a solid line 
we just cut the soft line to remind that no perturbative expression can be 
used for this line. Three relevant graphs are given on fig. 3 a,b,c. The full 
contribution of the 4-quark operators to the transition operator is 



Here 



— + + ^s- (22) 

Ld = 487r''ro^{Ki{cr^c){dr^d) + K2{cr^d)){dr^c)), (23) 

(24) 



A = 487r2ro^(ir3(cr^c - |c7^75c)(«r^xi) 
+i^4(cjr^Cfe - |cj7^75Cit)(Mfer^H')), 



where 



+K6{ciT^Ck - |Q7^75Cfe)(sitr^s*)), 



K, ^{cl + cl)/2, 
K, =(c^-d)/2, 

=-(c+ + c_)V8, 

= -(54 + c^-6c+c_)/8, 

=-(c+-c_)V8, 
Kq = -(54 + c2 +6c+c_)/8. 



(25) 



(26) 



The term Cd describes in the case of the meson the destructive inter- 
ference between the d quark produced in the charmed quark decay and the 
spectator d quark. The same operator in the case of describes the cd — > us 
scattering. The operators and Cg describe the interference between the 
spectator quarks in baryons and the -u or s quarks produced in the decay 
(this interference is destructive for the case of the u quark and constructive 
for the case of the s quark) . 

2.5 Hybrid anomalous dimensions 

Above the coefficients were presented for the operators normalized at rric. 
Generally speaking, these operators have the hybrid anomalous dimensions 
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[p!0| , 1^ reflecting the effect of the virtual momenta from rric to a typical 
hadronic scale /i. If we know the hadronic matrix elements of the operators 
normalized at rric we do not need to include the hybrid logs explicitly, they 
are already there. As we will see shortly this is the case with the operator 
^g{i^ = fTT'c) whose matrix element over mesons is expressible in terms of the 
measured mass differences while the matrix elements over and He vanish. 
The operators cc and cD^c have zero hybrid anomalous dimensions. 

As for the four-fermion operators, however, we have no independent in- 
formation on their matrix elements. Thus we have to use models - vacuum 
saturation in mesons and quark models in baryons. The question arises as 
to what choice of the normalization point would ensure the best possible 
validity of these models. 

As far as the strong interactions are concerned, is a completely foreign 
parameter. Therefore, it is absolutely natural to think that the models work 
best for a low normalization point, /i ~ a few hundred MeV. Then one 
should evolve the four-fermion operators down to fi, include explicitly the 
hybrid anomalous dimensions and only then apply the models. The effect of 
the hybrid logarithms is known in the leading-log approximation [^. 

Taking account of the hybrid logs changes the coefficients Ki introduced 
in eq. (|26|) as follows, 

^ K,r], (27) 

where 

K^vi ={cl + cl + l{l-^){cl-cl))/2, 
K2V2 = v^(4 - c2_)/2, 

K3V3 =-[(c+ + c_)2 + l(l-0^)(54 + c2 -6c+c_)]/8, 

7^47/4 = -v^(54 + d - 6c+c_)/8, ^ ^ 

= -[(c+ - c_)2 + i(l - v^)(54 + cl + 6c+c_)]/8, 
KgVq = -^{5cl + cl + 6c+c_)/8. 

Here 

The evolution from rric down to fi also adds new structures, for instance, 

AC, = 16n'—^To{cl - c'_)k'/'{k-'/' - l)(cVc)J^^ (29) 
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AC, = -8n'^K'/'iK~'/' - 1)(54 + c'_)icT,ec - ^C7^75t"c)j;. (30) 
Here 

The latter structures are similar to the penguin contributions in the weak 
charm decays. Since we neglected the penguins, it seems logical to discard 
the structures of the type ([29|), ( pOD as well. We refer the reader to ref. 
for a detailed discussion of the hybrid logarithms in the transition operator. 

Apart from the logarithmic renormalization of the currents and operators 
discussed above, hybrid logarithms also appear in the matrix elements of the 
4-quark operator when we use the factorization approximation to calculate 
them. We shall return to this point in the next section. 

Finally let us note that although the dimension 5 chromomagnetic oper- 
ator Og has non-zero hybrid anomalous dimension, we do not have to take 
into account the hybrid renormalization when we analyses its contribution. 
The reason is that the matrix element of Oq normalized at rric will be ex- 
pressed in the next section through physical quantities - the masses of the 
charmed hadrons . 

3 The hierarchy problem 

The most striking feature of the experimental data on the lifetime hierarchy 
obvious from the first glance at Table 1 is its wide spread, especially for 
the case of the baryons. This directly contradicts our intuition. Indeed, 
the known source of the differences in lifetimes between different baryons 
are the 4-quark operators. Although their coefficients are enhanced by a 
factor An^ their matrix elements are suppressed by a factor f'^/M'^ ~ 0.01. 
(Here fr, is the value of the decay constant of the D meson). Hence, one 
could expect the lifetime spread in the hierarchy to be at most of order 
one. Certainly one could not forsee that the lifetimes will span an order of 
magnitude, as experimental data show. The lifetimes differing by 50 to 100% 
was the prediction of the first works on this subject @, ^. In addition to this, 
the analysis [§], ^ had other problems. For instance, it was obtained that 
the lifetimes of A^ and are the same irrespective of the absolute value 
of < >. This result seems to contradict the experimental data. Some 
progress in understanding of the latter problem was achieved in ref. 0, where 
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it was noticed that the relative weights of different 4-quark operators are very 
sensitive to the Wilson coefficients c+ and c_. The use of the modern value 
of these coefficients corresponding to Aqcd ~ 250 — 300 MeV substantially 
improves the situation with Ac/H+ compared to the earlier results of ref. [^]. 
Still the problem of the wide spread remains. Can we explain it from first 
principles? The remainder of this talk will be devoted to an explanation (or, 
more precisely to speculations) why such a hierarchy can be expected from 
QCD (or, more precisely, how QCD must work so that the theory will be 
able to explain the experimental data). 

4 What Can Be Said about Mesonic and Bary- 
onic Matrix Elements 

Unfortunately, the issue of the matrix elements, the most underdeveloped 
element of the whole procedure at present, can not be deferred indefinitely. 
We have to address it, realizing, of course, that now we are not going to enjoy 
the same degree of theoretical control we had previously. 

For mesons the situation is not that bad, though, as will be seen below. 
At the same time in baryons we have to use much more shaky models. 

For each operator in the OPE we have discussed we consider first mesonic 
matrix elements and then, later on, the matrix elements over the baryon 
states. 



4.1 Mesons 

The first operator of dimension 5 is the operator caGc. The matrix ele- 
ments of this operator between the mesonic states were discussed in detail 
in ref. |2^. Remarkably, this matrix element can be determined in a model 
independent way, using the HQET symmetries: 

^ < D\ciaGc{fi = m,)\D >= /i^ = -{M\D*) - M\D)) ~ 0.4GeV^ 



(31) 

The next relevant matrix element is the operator of the kinetic energy 
^ < D\cD^c\D >= ~ 0.5GeV^ (32) 



2M 
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Here we used the QCD sum rule results pT |. 

We now move to the dimension 6 operators. Just one remark to demon- 
strate that Oe does not play a role. Indeed, assuming that the chromoelectric 
field is of the same order as the chromomagnetic field, we obtain that the 
matrix element of this operator is suppressed by a factor ~ /i^/m^ ~ 0.4 in 
comparison with the contribution of dimension 5 operators. The correspond- 
ing Wilson coefficient is 

ce = ^To/ml 

The contribution of this dimension 6 operator is less than 10% of the contri- 
bution of the dimension 5 operators and can be safely neglected. 

Let us turn now to the four-quark operators. The simplest approach 
allowing one to estimate the matrix elements of these operators is to use fac- 
torization. Consider first the operator cT ^^ddT fj_c. Note first that this matrix 
element contributes only to the D^, but not to the meson. This ma- 
trix element can be estimated using factorization, or the vacuum saturation 
method: 

< D+\cT^ddT^c\D+ >= flMl (33) 

where is the D meson decay constant, 

< Q\cT^d\D >= foip^. (34) 

The operator cT ^ddV and fr, are both normalized at one and the same 
point, rric. If we would like to proceed to fr, normalized at /i the correspond- 
ing hybrid logarithm should be inserted. 

We pause here to make a remark which may be very important numeri- 
cally. The question is which coupling fr, is to be substituted in eq. (|33|). The 
physical coupling fo contains in itself corrections that are known to be 
very significant. Therefore, if we want to have a systematic 1/mc expansion 
we must use a static value of /d, with the 1/mc terms removed. The differ- 
ence between the physical and the static value of fn clearly manifests itself 
in our expansion for the widths as an effect 0{m~'^). Were the expansion 
well convergent whatever value is substituted would not matter. For the ac- 
tual c quark this matters a lot, since the physical value of /d is ~ 170 MeV, 
while the static Fo ~ 400 MeV, i.e. ~2.5 times larger than the physical 
one [^, The difference between /I, and F|) is humongous. The static 
coupling constant is defined as follows 



< {)\dT^h\D >= iFDvf'JMD/2. (35) 
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We will use capital F to emphasize the distinction between the physical and 
the static constants. Fd defined by eq. ( p5|) depends on the mass of the 
heavy quark as 

Fd ~ mO)\VV^c (36) 

modulo hybrid logs. Here ip{0) is the value of the heavy meson wave function 
at the origin that does not depend on the mass of the heavy quark. If 
the value of matrix element ( p3|) calculated with Fd is substituted in our 
OPE we see that the corresponding contribution is pure compared to 

perturbation theory. 

On the other hand, the leptonic decay constant fn defined using eq. (P^) 
contains all orders of 1/ m^. Indeed, the relativistic Dirac field c{x) describing 
the c quark in QCD is connected with the the c quark field in HQET as 

c(x) = (37) 

Here D = ^^Da- Consequently, the static and physical decay constant are 
connected as 

/D = FB(l + Ci/m, + ...), (38) 

Ci is a constant independent of the heavy quark mass. It can be (and was) 
determined in the QCD sum rules ||2^ and on the lattices [^]. It turns out 
that the contribution of the term proportional to Ci in eq. ( |5BD is approx- 
imately 50% of the contribution of the leading term and has the negative 
sign. 

Only a systematic analysis of all dimension 7 operators can tell us what 
value - physical versus static - is more relevant to the problem. In the 
absence of such an analysis we will just speculate without even pretending 
on seriousness. 

Other four-quark matrix elements can be determined also using factor- 
ization and the Fierz identities. In particular, 

< D\cV^t''ddV^ec\D >^ 0. (39) 

For the product of two color-singlet brackets factorization becomes exact 
in the limit of the infinite number of colors. 

The natural question is how well factorization works in the real world. 
The corresponding estimate is completely analogous to the estimates of the 
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D — D mixing parameters where factorization works with the accuracy of 
order 20%. This means that for mesons we can determine matrix elements 
in a model independent way with a sufficiently high accuracy. 



4.2 Baryons 

We now move to baryons. Consider first the dimension 5 operators. We start 
from the operator caGc. The matrix elements of this operator between the 
members of the antitriplet of the charmed baryons are zero. Indeed, this 
matrix element is proportional to correlation between the spin of the heavy 
quark and that of the light cloud. The latter spin in the baryons - members 
of the antitriplet is zero. 

An estimate of this matrix element for the baryon belonging to the 
sextet of the charmed baryons (and so the light cloud there has spin 1) 
can be carried out in the same way as for mesons. Recall that a part of 
the Hamiltonian of the heavy quark effective theory responsible for the spin 
splittings between the hadrons has the form 

» hcBBhc , , . 

Since the chromomagnetic field is due to the light cloud, it is quite obvious 
that 

B = ^5^1.0., (41) 

where C is a constant and 1/2 is introduced for convenience. Then we easily 
obtain that 

AM2(Qf/ 2)) =< >= \c (42) 

and 

AM2(f]c) =< ^c\hc^BK\^^ >= -C. (43) 

Here n^?^'^^ is the corresponding member of the baryon sextet with the spin 
equal to 3/2. As a result, the matrix element of the operator Oq is given by 

— ^ < n.iOcl^c >= - < n,\h,aBhc\n, >= l{M^{n^^/^^) - Wisi,)). 

(44) 
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The mass difference M(fif/2)) _ M(fic) can be determined usine; nonrela- 
tivistic quark model ||25| and is ~ 50 MeV. 

Consider now the matrix element of the operator O-,^ = —cD'^c. The 
matrix element of this operator is nothing else than the average momentum 
square of light degrees of freedom. In ref. it is argued that the matrix 
elements of O.,^ are essentially independent of the particular hadronic state 
and are approximately the same for all c containing hadrons. In the absence 
of better ideas we accept this statement for orientation. If so, one can use 



the estimate of this matrix element from ref. |2T[ . 

The most difficult task in the baryon sector refers to the matrix elements 
of the four-quark operators. Unfortunately, in this case we do not have such 
a strong tool as factorization applicable in the case of mesons. In fact, at 
present we do not know any model-independent approaches permitting one 
to determine or even to estimate the matrix elements over the baryons. In 
order to say at least something in this case we have to turn to different models 
that are not based on first principles and are related to QCD only remotely, 
if at all. Consequently, the estimates of the matrix elements will be beyond 
theoretical control. 

In the literature analysis of the relevant matrix elements over the baryonic 
states was carried out only in the nonrelativistic quark model (NQM) and 
its close relatives, such as the bag model. The former will be discussed at 
length. Numerically the bag model predictions turn out to be rather close. 
As a matter of fact, essentially one feature of NQM is crucial - the absence 
of the spin-spin terms in the matrix elements of O^q over the triplet baryons 
(for brevity this clumsy combination of words is substituted by a shorthand 
'spin-spin interactions'). The overall normalization of the matrix elements 
does not matter much since effectively it is adjusted in an ad hoc way anyhow 
(cf. e.g. a 'static' prescription for M^^ — below). 

Let us recall some basic features of this model developed by De Rujula, 
Georgi and Glashow long ago. It is assumed that the hadron consists of 
2 (meson) or 3 (baryon) constituent quarks. The quarks in this model have 
'effective' masses, e.g. rric ~ 1.5 GeV, rriu ~ 0.34 GeV, etc. In order not to 
mix them with the current quark masses we denote them as m*. The quarks 
are described by nonrelativistic spinors and are bound by a nonrelativistic 



■^In NQM non-relativistic limit is assumed for both, c quark and the hght ones from 
the cloud. 
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potential modified by hyperfine interactions. The hamiltonian is equal to 

H = L{n,r2...) + SK + ^) + ^{aQiQj + ka,)S,j. (45) 

Here m* are the constituent quark masses, Pi are the quark momenta, and 
Tj are their coordinates; Qi are the quark charges, is the QCD coupling 
constant, and Sij denotes a relativistic hyperfine interaction. The constant 
k = —4/3 for mesons and k = —2/3 for baryons. The potential L is a 
universal quark-binding potential. The QCD interactions due to the gluon 
exchange are modeled by the Breit-Fermi interaction Sij assumed to have 
the same form as the electromagnetic Breit-Fermi interaction Sij We refer 
the reader to the original paper |2^ for a detailed discussion of this model 
and explicit (but rather complicated ) expressions for Sij. Using the mass 
dependence of the hyperfine interactions it is possible to show that in this 
model the absolute values of the squares of the wave functions of the baryon 
and the meson at the origin are connected as follows: [K 



l^(0)i?¥ ^ 2m:(M^.-M^.) 

mO)D\^ iMD-ml){Mn-Mn*)' ^ ^ 

Following refs. ^ ^ let us try to extract as much information as we can 
from the nonrelativistic quark model concerning the 4-quark operators under 
investigation. As we already noted above, the price will be the appearance in 
our equations of some parameters which are not connected in any way with 
the basic QCD, such as a constituent quark mass. Consider, for example, 
the Ac baryon. 

We begin from the matrix elements of the operator cT^cdT^d. These 
matrix elements are proportional to |V'(0)^^'^^P, the square of the baryon 
wave function at the origin. More exactly, this is the probability for c and d 
quarks to meet. The coordinate of the third quark is integrated over. From 
now on the superscript cd will be omitted. The relation (|46|) connects the 
squares of the wave functions in the baryons and mesons [^, ^ and will be 
systematically applied to make estimates in baryons. Note that in the limit 
TTT-c — *■ cxo we are interested in 2M£)(M£). — M^,) = /i^ x 4/3 and Mj.+ — 
tends to a constant value. Consequently, 

3(M^+-M. + )m* ^ o 1 o ./n 

\m^f- 2 ' \mD\\ m)D? = -FiMnK-"\ (4?) 
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Recall now ^ that the matrix element of the operator cTf^cdTf^d over 
the baryon He is equal to 

^ < H,\cT^cdT^d\H, >= \^hM\\^ - 45,5^,) (48) 



where S denotes the spin operator. In baryons from the antitriplet - Sc and 
Ac - there is no correlation between the spin of the heavy quark and the 
spins of the light quarks. Hence the spin part of eq. (^) is zero. Using the 
expression for |'?/;(0)£)| in terms of F]j and M]j we immediately obtain 

< A,|cr^cJr^d|A, >~ - M^+)m:F^M^«:-^/^ (49) 

In the same way we find for the sextet Qc baryon 

1-5 
— — < riclcT ^cdT ^d\ilc >~ — 
2Mn^ 6/iG 



< Qc\cT^cdT^d\nc >~ ^{Mj.+ - M^+)mlFlMDK-^'^. (50) 



Here we assumed that the spin splittings like Mj.+ — do not depend on 
flavor and are the same for all charmed baryons. 

The matrix elements of other operators over A^ can be expressed via 
matrix element (^) using the Fierz identities and antisymmetry in color of 
the baryon wave functions. In particular, using the color antisymmetry of 
the baryon wave functions we obtain that || 

< Ac\ciT^c^djTf,d'\Ac >=< Ac\cT ^ddV f,c\Ac >= - < A^l cF^c Jr^d|Ac > . 

(51) 

Similar relations can be found for all other charmed baryons. 

Note that if we want to pursue a systematic expansion in l/rric we must 
use the static value of the meson decay constant in all expressions above. 

5 Static Versus Actual fq and Speculations 
on Higher Order Corrections 

Above we have carried out a systematic expansion of the inclusive width 
in l/mq including the contributions of all operators with dimensions up to 
6. Let us consider our estimates of dimension 6 four-quark operators more 
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carefully. As it was already mentioned if we want to have a systematic 
l/mg expansion we must throw away the part of fn that contains the terms 
suppressed by powers l/mg and use the static meson decay constant in 
all our estimates. If we want to deviate from this procedure we must take 
into account the terms of dimension 7 and higher. 

A natural question is whether we have to deviate from this procedure. 
After aU, it is reasonable to expect that the contribution suppressed as 
will be unimportant. Unfortunately, this is not the case. If we combine the 
formulae of section 2.5 with the estimates of the 4-quark operators obtained 
using Fd, we obtain that they give the contribution to the total width of 
meson which is negative and ~ 3 x the width of meson. This is clearly 
nonsense, meaning that the higher dimensional operators play an important 
role in the hadronic width of D'^ mesons. Can we model in any way the 
contribution of the higher-dimension operators? Can we at least argue what 
they must show up in order to have sensible results? 

First, consider the mesons. We shall argue that in the meson case there 
is at least a set of higher-dimension operators whose effect is to renormalize 
Fd to /d. 

Indeed, consider the operators appearing in the expansion of the four- 
quark operators in terms of the operators of HQET. Recall that heavy quark 
currents like cTc have the following expansion : 

(2MD)~^crc^ hcThc 

. t: (52) 

+{i/2MD)ihcrDhc -hcD The) + 0{l/ml). 

Here I) = ■j°'Da. The four-quark operators we are interested in also have a 
similar expansion, 

^cT^cqV^q = {hcT f,hcqT ^.q 

- T (53) 
+{i/2MD){cr^Dhc -hcD rhc))qr^q + 0{l/ml). 

The matrix element of the operator (53) is naturally expanded in terms of 
the HQET operators. The matrix element of the first term in this expansion 
is equal to 

< D\hcT^KqT^q\D >= F|,Mb/6. (54) 

We used the Fierz transformation and factorization to estimate the latter 
matrix element. 
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By the same token, assuming that factorization still holds for the next 
terms we obtain that the matrix element of the second term in the expansion 
( |55D is such that 

-J- < D\cT^cqT^q\D >= \mdFI{1 + 2Ci/Md + ...)• (55) 

The right hand side of the latter equation is nothing else than the first two 
terms in the l/rUc expansion of f^M^. 

The above arguments show that it is natural to assume that within fac- 
torization the role of the higher dimensional terms is to renormalize to 
Jd in calculating the 4-quark matrix elements. Unfortunately, the situation 
is not that simple. First, there are other terms of higher dimensions beyond 
the ones that appear in the HQET representation of the four-quark opera- 
tor. Second, even among the operators that appear in this representation 
there are " nonfactorizable" ones, that can not be reduced to the form of the 
product of the light quark current - heavy quark current. We can only hope 
that future investigation will clear the situation. 

Note that the need of renormalization of Fd directly follows from the 
hope that our theory has sense in the charmed family. Indeed, if we use 
in the estimates of the 4-quark operators in mesons, we obtain that the width 
of meson is negative (~ —2 widths of the meson), a clear nonsense. 
On the other hand, if we use fo, we obtain something more reasonable and 
in accord with the data. 

Following the above discussion we conjecture that the effect of the higher 
dimension terms is to renormalize the leptonic decay constant of the heavy 
hadron in the estimate of the 4-quark matrix elements from its static to its 
real value provided that factorization is applicable. We use the 'real world' 
value of the leptonic decay constant while doing calculations in mesons. 

You can view the speculation above (/^ versus Fd) as just an attempt 
to disentangle the overall normalization of the mesonic matrix elements from 
that of the baryonic matrix elements, a connection inherent to NQM, if you 
are not inclined to take it more seriously. 

Consider now the case of the baryons. In this case we can not repeat 
the above calculations - for baryons factorization is not applicable. We can 
not estimate the contribution of the operators appearing in the second line of 
(53). We, thus, have to be satisfied with three first orders in l/mg expansion, 
and have no alternative but to use the static constant Fd in our calculations 
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assuming that here the higher order terms do not play as drastic role as in 
the meson lifetimes. Quite remarkably, we will see later on that then we do 
get the experimental pattern of the lifetimes. 

We end this section by reiterating that the preasymptotic power terms 
must be much larger in baryon lifetimes than in meson one, and this can 
be achieved by using for baryons and fo for mesons. The difference is 
formally of higher order in l/rric - it is necessary to have enhanced terms 
(!?(m~^) in mesons and suppressed in baryons. 

An obvious difference between the baryon antitriplet and mesons is the 
presence of the strong spin terms in the latter. The same strong spin terms 
are present in the Qc baryon. It is, perhaps, possible to conjecture that an un- 
known mechanism that distinguishes between the mesons and the antitriplet 
baryons is somehow connected with the spin terms. If so, a natural conse- 
quence might be that the charmed heparin Qc has the lifetime around that 
of D^, i.e. larger than the lifetimes of or It will be very interesting to 
determine this lifetime experimentally. 



6 Numerical Analysis in the Charm Family 

We move now to numericals. We shall consider first the case of the charm 
family. The following numerical values of different parameters are accepted. 
We use the values c+ ~ 0.74, c_ ~ 1.82 corresponding to the currently 
accepted value of Aqcd ~ 300. We use fc ~ 170 MeV for the 'real world' 
value of fr), and Fd ~ 400 MeV for the static leptonic decay constant. 



both normalized at p^ , [23|] . We choose the normalization point ~ 0.5 
GeV for the estimates of matrix elements of the four-quark operators. The 
corresponding value for k is k ~ 3.1. We use the numerical values rric ~ 1.4 
GeV, ~ 0.4 GeV2, and ~ 0.5 GeV^. 

We also use Mj.+ — Ma^ = 400 MeV, not the real world value M-^+ — 
= 200 MeV, in the systematic l/m^ expansion, since, as it follows from 
QCD sum rules p3, this mass difference goes to ~ 400 MeV in the limit 



of the infinite mass of the heavy quark. (Remember, this is just an overall 
normalization of the NQM matrix elements, see eq. (46)). For the constituent 
quark mass we use m* ~ 0.35 GeV. 

In the naive parton model (the dimension 3 operator Oq) all hadrons have 
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the same width 

r = N,7]To ~ 4.5ro. (56) 
They also have the same semileptonic widths 

Ts/i ~ 2ro. (57) 

The total width in this model is ~ G-SFq. 

Consider next the contribution of dimension 5 operators. We begin from 
the charmed mesons. These operators are spectator blind and contribute in 
the same way to and D~^. Using the explicit formulae for the matrix 
elements of Oq, Oq, C'tt and the Wilson coefficients Co, Cg, C-,, it is easy to 
obtain that the dimension 5 terms from Oq, Oq and 0-,^ shift hadronic widths 
of D mesons by 

Ar(^) ~ 2.5ro. (58) 

Let us now consider in more detail the case of meson. Recall that its 
width also gets contribution from annihilation mechanism ~ O.OFo |^|, and 



that its semileptonic width is suppressed by a factor ~ 1/2 |1|, |T9[ . Then the 
hadronic width of meson is 

ThMD') ~ 8To. (59) 

and its total width is 

r(D°) ~ 9To. (60) 

Of course, all these estimates have only qualitative character although we 
certainly are always glad when we see the right tendency. To give a feeling of 
the uncertainty, note that if we use ~ 1.35 GeV, the width of increases 
by ~ 10%. 

This is all with the D° meson, since the 4-quark operators of dimension 
6 do not contribute to its width. 

Let us parenthetically note that Dg is in the same situation as and is 
predicted to have the same width in the limit when we neglect SU (3) breaking 
effects. A tiny difference between Dg and is exhaustively discussed in a 
recent work ^9], and we will not dwell on this point here. 

For the baryons (except flc) the contribution of dimension 5 operators is 
due only to operator and is ~ — 0.5ro. The chromomagnetic operator 
Og proportional to the correlations between the spin of the heavy quark and 
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the spin of the hght diquark does not show up since the spin of the latter is 
zero for the baryon antitriplet. For Qc we estimated the contribution of this 
operator in section 4.2. 

We now move to dimension 6 operators. The values of the Wilson coeffi- 



cients KiTji, eq. (|28|) , are: 



Kirii ~ 2.28, 
~ -2.45, 

K3V3 ~-0.88, , . 

Ki7]i -0.45, ^ ' 

-0.31, 
K^riQ ~ —2.65. 



Consider first the charmed mesons. The contribution of the 4-quark op- 
erators in the width of meson is zero. For meson the contribution of 
the 4-quark operators is given by 

rW = 48.=r„(:^ + /f.,.)Aj_. (62) 

Here takes into account the hybrid logarithms that appear in /d when 

we move between rric and a low normalization point. Substituting the numer- 
ical values of the parameters given above and using the "real world" value of 
the constant /d ~ 170 MeV we obtain that the contribution of the dimension 
6 operators into the width of the meson is 

Ar(^)(D+) ~ -4ro. (63) 

For the total and hadronic widths of the meson we immediately obtain 
(taking into account the suppression of the semileptonic width mentioned 
above, and the contributions of Oq and the dimension 5 operators that are 
the same as for meson) 

rhadr(/^+) ~ 4ro, r(D+) ~ 5ro, (64) 
r(/^°) . ^ 

This result in an agreement with the experimental pattern, as can be easily 
seen from Table 1. 
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This agreement, however, must not be taken too hterally. Indeed, in 
the latter estimate we neglected the difference between Mo and rric in the 
denominator of eq. (|52D. Were we distinguishing between M^) and we 
would obtain, instead of the factor /|,/M|, in eq. (|62D, the factor f'^jMo/w?^ 
that will give a huge contribution of dimension 6 operators, close to — TFq, 
reducing the width of almost to zero. This once again shows the im- 
portance of the contribution of higher- dimension corrections. In our case it 
seems logical to substitute rric by Mr, in the estimates since we use the real 
world value of fo- If we used the systematic l/rric expansion and the static 
value of leptonic decay constant Fr, we would obtain an unphysical negative 
width of meson. Alas, the charmed hadrons are too light! 

Let us now consider the baryons. Using the explicit expressions for the 
baryon matrix elements in the NQM it is straightforward to get the following 
expressions for the matrix elements of the dimension 6 operators: 

Wu =j^<H,\C^\H,>={K3m-K^7]^)6TB, (66) 

Here we neglected the penguin- type structures. The width 6Tb is determined 
by the following expression 

o m*,F'lMD{M^+ - Ma ) 

The value of STb depends on whether we use Fd and what value of {M-^+ — 
M\J is substituted. If the static mass difference, 400 MeV, is substituted, 
along with the static Fd, we immediately obtain that the contributions of 
the different 4-quark operators are: 

~28ro, 

=~-8ro, (68) 
~22ro. 

For the contributions of the 4-quark operators into the hadronic widths of 
the charmed baryons we find 

r£L(A+) = w^, + w^. ~ 2oro, 

rSUs^) =W^ + Wsr^UTo, (69) 



27 



Let us now consider the full baryonic and mesonic widths. Taking into ac- 
count that the semileptonic widths of the triplet baryons are approximately 
the same as in the naive parton model (no Oq contribution) and assembling 
together operators Oq, and 4-quark contributions we immediately obtain 

r(A+) = 26ro ~ 5r(D+), 

r(s+) = 2oro ~ 4r(D+) ~ 2r{D^), 

r(SO) = 56ro ~ 11T{D+). 

These results are in a qualitative agreement with the experimental pattern 
in Table 1. 

If instead of the static Fjj we used /d, we would get numbers completely 
contradicting the experimental pattern and leading to too narrow a span of 
the lifetimes. 

We conclude that the scale of the matrix elements of the four-quark op- 
erators must be enhanced in baryons compared to mesons. Within NQM 
we achieve this enhancement by applying a prescription Jd Fd- From 
theoretical point of view the obvious difference between the mesons and the 
antitriplet baryons is the absence of the interactions between heavy quark 
spin and that of the light cloud in the latter. Is this the underlying reason 
for the different behavior of the heavy quark expansion in these two cases? 
There exists one charmed baryon, namely f2c, that has a strong spin-spin 
interaction like in the mesons. If we use the "static" prescription, as for an- 
titriplet baryons, we obtain that the width of this baryon is ~ 16r(D+). If 
we follow the meson example and use the prescription, we obtain that its 
width is ~ 2V{D^). Depending on the prescription used VLc can be both the 
most short-living and the most long-living among charmed baryons. There 
is yet no experimental data on the width of this baryon. Such data can be 
crucial to establishing the QCD mechanism of the hierarchy formation. 

Note the main distinctions between our analysis and the analysis carried 
in refs. |^, 0. In comparison with the first two works, we used different 
values of the Wilson coefficients. The ones used in correspond to Aqcd ~ 
100 MeV. New values of the Wilson coefficients change in the right way 
the relative contributions of different four-quark mechanisms (proportional 
to Wd,Wu,Ws, respectively) in the inclusive widths of the baryons. We 
also discarded a rule of discarding accepted in Instead, we took 

into account the contribution of dimension 5 operators. Next, following the 



authors of ref. HI we used the relation E6| between the absolute values of the 
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wave functions of baryons and mesons. Finally, we used the static value Fn 
for the estimates of the baryonic matrix elements. 

In comparison with the analysis we included the hybrid logarithms 
and the contribution of dimension 5 operators. 

7 Estimates for Beautiful Hadrons 

The analysis for beautiful hadrons proceeds in exactly the same way as for 
charmed ones. Moreover, we expect our predictions to be more accurate, 
since the expansion parameter is ~ 0.1. 

The only difference with the charm case is the presence of the phase space 
factors: /o = 1 — 8r + 8r^ _ _ 12r^lnr for the cc operator and (1 — r)^ 
and (1 — rY for Oq and O^, operators respectively. We use the 6-quark mass 
mft ~ 4.8 GeV and the values /s ~ 150 MeV, Fb ~ 250 MeV fS^, 0. The 
normalization point /i ~ 0.5 GeV, corresponding to k ~ 5.42 is chosen. The 
values of the Wilson coefficients are 

c+ = 0.84, c_ = 1.41. 

The question of the lifetime hierarchy is especially interesting, since in 
this case we expect the inverse heavy quark mass expansion to work much 
better than for the charmed hadrons. All calculations become more reliable. 
It will be of great interest to see whether the approach becomes quantitative 
for the beautiful hadrons. Recall that we speculated that fo is more relevant 
for mesons while -Fp is more relevant for baryons. The difference between 

and F| is still quite noticeable - about 40%. The same regularity is 
expected to show up in the beautiful hadrons. Unfortunately, at this point 
we do not have any reliable experimental data, so we just consider theoretical 
predictions that follow from the use of the same rule for beauty as for charm. 

Using the values of the parameters given in the beginning of this sec- 
tion, we immediately obtain rather large lifetime differences between differ- 
ent beauty mesons. In particular, we obtain that the lifetimes of charged and 
neutral B mesons can differ by as much as 3 to 5%. The width of the charged 
meson T{B~) decreases due to the interference effects. On the other hand 
we obtain that the width of A;, can be as much as 10% bigger than the width 
of the B°-meson and is approximately equal to the width of the baryon. 
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The width of H° baryon is ~ 3% smaller than the width of meson. It will 
be extremely interesting to check these predictions experimentally. 

Note that the significant difference with the previous estimate of the 
beauty lifetimes comes because we use the value of /b ~ 1-5 larger than 
Voloshin and Shifman, as it follows from the recent QCD sum rule analysis 

8 Conclusions 

In this talk we addressed the question whether it is possible to 'explain' the 
observed hierarchy of the hadronic widths in the charmed system in QCD 
on the basis of the heavy quark expansion. Strictly speaking the charm 
family seems to lie below or, with luck, at the border of the domain where 
the expansion may be useful, which makes the task absolutely challenging. 
We speculated that there is a significant difference in the structure of the 
preasymptotic expansion in the meson and baryon systems. Terms of higher 
dimensions show up in a completely different way in mesons and baryons. 
This situation can be effectively modeled by taking into account the first 
orders of the preasymptotic expansion (up to C(m~^)) and by using the real 
world value of the constant fo in mesons and the static in baryons. If 
such a picture really exists the problem of a wide spread of the total widths 
in the charmed family is solved, at least at the qualitative level. 

Unfortunately, the origin of this basic difference between mesons and 
baryons is still unclear. Since the obvious difference between mesons and the 
antitriplet baryons is the absence of the spin-spin interaction in the latter 
(and its presence in the former), it will be extremely important to get a 
reliable data on the lifetime of the sextet baryon Qc- This is the only baryon 
where we expect that the spin-spin interactions play a major role, as in 
mesons. 

Assuming that the same mechanism works also for beauty hadrons, we 
listed in the last section some estimates for beauty hadrons. It is extremely 
interesting to check them in current experiments. 

Finally we must stress that our estimates are just provocative. Our task 
here was , first of all, to formulate the problem of the hierarchy, to see what 
is known in the current literature about it, and whether it can be solved 
using available methods. We have speculated how QCD must work in order 
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to explain the hierarchy. The next question is whether it really works that 
way. 
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Experimental Data on 



Table 1 

Weak decays of Charmed Hadrons 



hadron 


lifetime 


D+ 


~ 10.6 X 10-^3 s 




~ 4.2 X 10-^3 g 


Ac 


~ 2 X 10-^^ s 




- 3 X 10-13 g 




~ 0.8 X 10-13 s 
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Figure Captions. 

Fig. 1: The imaginary part of this graph determines the hadronic width 
of all c containing hadrons in the parton model approximation. The closed 
circle denotes the effective weak lagrangian. 

Fig. 2: The influence of the soft modes. A soft gluon (absorbed by the 
light cloud of D) is marked by x. 

Fig. 3 a,b,c: The diagrams giving rise to O^q in T. The crosses mark 
soft d, u or s quark lines 
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